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ИНФОРМАЦИОННОЕ МОДЕЛИРОВАНИЕ р 
СИСТЕМНОЙ ДИНАМИКИ ПОРШНЕВЫХ ДВИГАТЕЛЕЙ ВНУТРЕННЕГО СГОРАНИЯ 


Детализирована модель рабочег О процесса поршневого двигателя. Выполнено сопряжение динамики внут- 
рицилиндровых процессов с внешними условиями — нагружения, теплообмена и топливоподачи. Разработан 
дополнительный модуль, который позволяет удовлетворительно описывать кинетику окислов азота. Частич- 
но решена задача расчета химического состава отработавших газов в различных режимах, включая пере- 
ходные. По сравнению с прототипами модель обладает повышенным химическим разрешением, что позволя- 
ет с ее помощью уточнять теоретические оценки содержания отдельных токсичных веществ в отработавших 
газах. 
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Введение. Поршневой двигатель внутреннего сгорания (ДВС) является основным агрегатом со- 
временных автотранспортных средств, поэтому топливная экономичность и токсичность отрабо- 
тавших газов (ОГ) имеют первостепенное значение. Отсюда растущий интерес к информацион- 
ным моделям рабочего процесса ДВС как инструменту интерпретации данных, оптимизации кон- 
струкции, средству поиска альтернативного моторного топлива [1, 2]. 

При моделировании внутрицилиндровых процессов приходится преодолевать «проклятие 
размерности», проблему разнородности и разномасштабности процессов, аккуратно описывать 
«обостренные» режимы и формирование структур. Важно отработать методики агрегирования 
фрагментных моделей таким образом, чтобы результат «склеивания» наследовал регуляризую- 
щие свойства ключевых элементов [3]. 

Попытки рассчитать токсичность поршневого ДВС при помощи компьютерных моделей 
[4-5] лишь отчасти оправдали ожидания авторов: хорошо воспроизводя выбросы СО, углеводо- 
родов С,Н,, сажи и сернистых соединений, они давали на два порядка заниженные концентрации 
МО в ОГ. Поскольку в первых моделях внутрицилиндрового горения в двигателе азотная кинетика 
не учитывалась, проблема рассчитать содержание МО в выхлопе связывалась с химической не- 
равновесностью рабочего тела, которое полагалось параметрически однородным. Впоследствии 
кинетическая схема [5] была дополнена реакциями с участием азота, но для расчета концентра- 
ции МО в ОГ этого оказалось недостаточно. 

Ключевые идеи, допущения, математическая формулировка задачи. 

Отработана следующая гипотеза. Если макроскопически трехзонные кинетические модели не по- 
зволяют адекватно воспроизвести токсичность ДВС по МО, причиной может служить только про- 
странственное переосреднение, т.е. замена истинного поля параметров Х(х, #) средними по про- 
странству величинами, например, <Х(#)> или Х(<у>, 8. Это приводит к тому, что важные для це- 
лей моделирования связи не отражаются моделью. В математической формулировке задачи ис- 
пользованы обозначения: Х= {ст, со, ..., Г} — псевдовектор, объединяющий концентрации компо- 
нентов с;, см, и температуру Т, К; г = #1, г», гз} - пространственная вектор-переменная, 
{м, м, м}; Е- время, с. 

Информационная модель (ИМ) упрощается, если в системе имеются и выявлены малые 
или, наоборот, большие параметры [3]. 

Схематизация горения в двигателе: модель «реакции-перенос». Основополагающие 
предположения усовершенствованной модели унаследованы от прототипов [4, 5] и формулиру- 
ются следующим образом: 

— топливно-воздушный заряд химически и термодинамически однороден; 

— впрыскиваемое топливо монофракционно, а его условная химическая формула С!.Нэо 
(дизельное) или СзНаз (бензин); 

— химические аспекты горения хорошо укладываются в кинетическую схему, ограничен- 
ную элементарными реакциями между С-Н-О-содержащими веществами с тремя и менее атомами 
углерода, а также несколькими брутто-реакциями крекинга исходных углеводородов; 


509 


Технические науки 








— турбулентные эффекты сводятся к ограничению скорости элементарных химических 
процессов темпом переноса реагентов в зону тепла из зоны пламени; 

— испаряется жидкое диспергированное топливо диффузионно; 

— теплообмен излучением незначителен. 

Поскольку диффузия, теплопроводность и реакции — равнозначные «игроки» в процессе 
формирования пламенных структур, кинетические члены в уравнениях непрерывности и энерге- 
тического баланса подверглись «мелкомасштабному» осреднению. Для учета макроскопической 
параметрической неоднородности топливно-воздушной смеси осуществлено разбиение заряда на 
три взаимодействующих зоны, внутри каждой из них рабочее тело полагается однородным. Ско- 
рость межзонного обмена термохимическими характеристиками определяется как «качество сме- 
сеобразования» с помощью полуэмпирической зависимости. 

Энергетическая открытость модели обеспечивается заданием реалистичных условий теп- 
лоотдачи со стороны заряда (горящей топливно-воздушной смеси) в стенки камеры сгорания 
и внешнего силового (или, точнее, моментного) нагружения ДВС. Теплообмен задается уравнени- 
ем Ньютона-Рихмана, а нагрузка — известными уравнениями кинематики кривошипно-шатунного 
механизма [6] под действием внешних сил, препятствующих свободному вращению коленвала. 
В модели эти силы имитирует линейная комбинация сухого трения и гидродинамического сопро- 
тивления движению, которое выбрано пропорциональным частоте вращения коленвала. 

В разработанной авторами информационной модели (ИМ) рабочего процесса ДВС химизм 
горения углеводородного топлива описывается совокупностью более чем 400 микроскопически 
обратимых элементарных реакций, кинетические характеристики которых взаимно согласованы 
в рамках общего термохимического потенциала реагирующей газофазной среды [7, 8]. Эта важ- 
ная особенность представленной ИМ качественно отличает ее от известных аналогов и прототи- 
пов, в динамике которых не исключены артефакты. 

Собственно внутрицилиндровый химизм, его связь с динамикой и теплообменом описыва- 
ется системой дифференциальных уравнений 

+800, , (1) 

где /- скорость изменения термохимических параметров заряда Х, {см?. с!, см? .с*",...,К.с\}; 

5- темп обмена термохимическими параметрами между зарядом и «резервуаром», {см?. с', 
Е НИ 

Система (1) подразумевает пространственное осреднение динамических переменных, при- 
чем в самом общем виде. Уравнения внутрицилиндровой динамики ДВС и сопряженных процес- 
сов, осреднением которых получены (1), подробно описаны в работах [7, 8], а основная идея ме- 
тода их осреднения описана ниже. Ключевой особенностью описанной модели является двухмас- 
штабное осреднение динамических переменных, объединяющее параметрическое осреднение ки- 
нетических уравнений на уровне пламенных диффузионно-реакционных структур (с учетом тур- 
булентного разлохмачивания последних) и осреднение на уровне квазиооднородных макрообъе- 
мов горящей топливно-воздушной смеси (ТВС) заряда. 

Микромасштабное осреднение эволюционных уравнений основано на следующих научных 
результатах. Распространяющееся по однородной горючей смеси пламя служит классическим 
примером автосолитонов [8]. Выделяющаяся в процессе горения среды внутренняя энергия «со- 
бирается» вследствие конкуренции нелинейности и диффузионного переноса в бегущую локали- 
зованную структуру (волну). Скорость распространения волны горения — среднегеометрическая 
величина между обратным временем химической реакции и эффективным коэффициентом диф- 
фузии вещества и/или тепла в горючей среде. Поэтому явной формой эволюционных уравнений 
динамики поршневого (одноцилиндрового) ДВС, осредненных по структуре/масштабу пламенной 
зоны, являются: 
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ЕЕ ЖЕРТЕТЬ 


п = та, п; = с; /р[молекул/г] - концентрация компонентов (р — плотность газа, г/см3); 
Ве - число Рейнольдса; Р;— внешние источники газообразного вещества (при впрыске); 
- доля теплового потока, не связанная с химическими реакциями; ти Г,- соответственно 
масса жидкой фазы и теплота испарения; и; — молекулярный вес /-го газового компонента; 
5„ и [— мгновенное значение площади поверхности теплообмена и характерного линейного 
размера камеры сгорания (КС); 7 и = - отношение радиуса кривошипа к длине шатуна и 
геометрическая степень сжатия; а, Т„, — усредненные коэффициент теплоотдачи, темпера- 
тура поверхности КС; О; Л = Л(С) и С(1) — коэффициент диффузии паров топлива в вакуум, 
средний диаметр капель при распылении форсункой и массовая скорость впрыска топлива; 
бу — символ Кронекера; К - универсальная газовая постоянная; / — приведенный момент 
инерции нагруженного ДВС; Р(ф), 5, и ’(ф) — соответственно мгновенное давление в камере 
сгорания, площадь поршня и имеющая размерность длины функция передачи поступатель- 
ного усилия во вращение (это слагаемое в уравнении (11) соответствует одноцилиндровому 
ДВС); ‹ф - угол поворота, град. ПКВ и О=аф/@- частота вращения коленвала; 
Е и ^ - коэффициенты сопротивления вращению; А — постоянная Больцмана; прочие обо- 
значения соответствуют принятым в теплотехнике. 

Фигурирующий в уравнениях (2)-(12) комплекс Г” агрегирует два параллельных канала 
химических превращений: / “"” — медленные реакции, протекающие синхронно во всех точках 
камеры сгорания как в реакторе идеального смешения, и /“ — быстрые процессы во фронте пла- 
мени, темп которых ограничивается турбулентной диффузией. Раскрывающие смысл (2)—(12) ал- 
гебраические соотношения уравнивают роли реакционного и транспортного каналов в совокуп- 
ном процессе образования/гибели 1-го вещества. В соответствии с (2)-(12) при идеальном пере- 
мешивании все вещества образуются равномерно в объеме камеры сгорания, а при сильно за- 
трудненном -— во фронте горения. Мерой перехода того или другого канала образования вещества 
в ранг доминантного служит степень турбулизации заряда. Предложенная техника позволяет дос- 
таточно аккуратно учитывать роль диффузионных процессов при горении в ДВС, используя обык- 
новенные дифференциальные уравнения в качестве базового математического инструмента. Ес- 





| 
РВе" и, | Х(и,р,Т+Т.,/2)|; 
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тественно, в пламени реакции протекают при повышенной температуре: ее скачок ограничен 
адиабатической температурой сгорания ТВС в локальных условиях ТГ; и в нашей модели полага- 
ется равным 7.:/2. 

Кроме того, поскольку нет универсальной связи между коэффициентом турбулентной 
диффузии Ди и числом Ве, в нашей модели использована зависимость О» — Ве”, промежуточ- 
ная между известными предельными приближениями, мелко- и крупномасштабных пульсаций 
Ри — Веи Ри 85 Ве’. 

Определенным недостатком предложенных здесь уравнений горения является отказ на 
данном этапе от осреднения скоростей реагирования по Фавру на масштабах пламенной зоны [9]. 
Оправданием такого упрощения нашей модели служит достаточно интенсивная турбулизация, 
следовательно, заметное размытие реакционной зоны в пространстве. Согласно специально вы- 
полненным оценкам, в условиях ДВС чувствительность констант скоростей реакций к изменениям 
температуры заметно превышает их чувствительность к плотности реагирующего газа. 

Если ограничить осреднение уравнений химизма только мезотранспортными процессами 
в пламени, концентрация ряда токсичных продуктов сгорания оказывается заниженной. Это каса- 
ется как продуктов неполного сгорания СО и С.Н, так и побочных продуктов окисления МО,. При- 
чина кроется в том, что выход этих продуктов определяется крупномасштабной параметрической 
неоднородностью ТВС, а не кинетикой гомогенных реакций. 

Для приведения модели в соответствие реальному положению вещей в ней использовано 
предположение о существовании в пределах заряда ТВС трех параметрически однородных зон: 

— с номинальным коэффициентом избытка воздуха ои=о,, в которой сгорает х-я часть топ- 
лива (0 <х < 1); 

— коэффициентом избытка воздуха а» = 0,7, где сгорает 1/2(1 — х) доля топлива; 

— коэффициентом избытка воздуха оз = 2, — 0,7, где сгорает остаток топлива. 

Итоговую концентрацию /го компонента находим при помощи соотношения 

<сР= и (1-х с + хе" и 1-х с, (13) 
где верхний индекс при с; соответствует локальному коэффициенту избытка воздуха. 

Для расчета качества смесеобразования х нами предложена формула 

х = Ш(телет/ ва) = Ш[Оиь п! (РУ Ата, (14) 
ГДЕ телем И Та — характерное время химического окисления и механического перемешивания заряда. 

Внесение коррективы в исходную модель значительно улучшило согласие рассчитываемых 
величин с;с известными из эксперимента тестовыми данными. 

Параметры модели 5х, / и Ве (соответственно Д‚шь), безусловно, являются конструктивно 
определенными функциями времени (. 

Развиваемая в соответствии со схемой (рис.1), модель позволяет не только надежно оп- 
ределять энергетические (интегральные и мгновенные) характеристики рабочего цикла различ- 
ных типов ДВС, но также содержание в ОГ СО, С,Н„, МО, ВСОН, воспроизводить тонкости образо- 
вания токсичных веществ на переходных режимах. 
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Энергетика + динамика + суммарная токсичность отработавших газов 


Рис.1. Эволюция модели рабочего процесса поршневого ДВС [4—8] 
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Авторами также ведутся работы по расширению номенклатуры токсичных компонентов 
в выхлопных газах за счет сажи и бенз(а)пирена С»оН:2. 
Результаты вычислительного эксперимента и их интерпретация. Адекватность созданной 
модели оценим, сопоставляя рассчитанные с ее помощью характеристики рабочего процесса ДВС 
и данные технического эксперимента. Соответствующие результаты для распространенного теп- 
ловозного дизеля 1Д49 приведены в табл.1. 
Таблица 1 
Некоторые интегральные, экстремальные и мгновенные характеристики 
рабочего процесса в дизеле 1Д49 в режиме номинальной мощности 


| Максимальный темп нарастания давления, атм/град ПКВ 3 
| Максимальная температура, К | 16300 | 
| Максимальная скорость роста температуры, К/град ПКВ | 610 | 
[Доля выделившегося тепла заряда до Рив», % | 490 | 





Эти данные свидетельствуют, что доступные для экспериментальных измерений парамет- 
ры рабочего цикла надежно рассчитываются на основании предложенного авторами кинетическо- 
го подхода. Кроме того, разработанная схема вычислений позволяет определять ряд не измеряе- 
мых, но важных для практики эксплуатации ДВС мгновенных и интегральных параметров рабоче- 
го процесса (выделены в табл.1 жирным шрифтом). 

Для целей оптимизации впрыска и газообмена в поршневом ДВС важны интегральные 
и мгновенные энергосиловые показатели рабочего процесса (давление и температура в камере 
сгорания), а также токсические характеристики заряда и ОГ. Разработанная для управления элек- 
тромагнитной системой клапана математическая модель внутрицилиндровых процессов позволяет 
рассчитывать определяющие параметры рабочего тела в динамике с большей (механические) или 
меньшей (токсические по оксиду азота) степенью достоверности. В табл.2—4 приведены результа- 
ты энергосиловых расчетов для трех представительных режимов двигателя ВАЗ-2108: скорость 
вращения коленвала составляет на этих режимах 1000, 2000 и 4000 мин`!, коэффициент избытка 
воздуха си = 1. 























Таблица 2 
Параметры рабочего процесса двигателя ВАЗ-2108 
Коэффициент избытка воздуха, 1 0,7 1,3 Среднее 
доля объема 0,8 0,2 0,2 значение 
1 2 3 4 5 
5,93 5,91 1,34 5,47 
Среднеиндикаторное давление, атм 6,47 5.82 1,31 5,89 
5,86 5.67 1.59 5,41 
16 15,9 3.6 14,8 
Индикаторная мощность, кВт 34 31,4 7,1 31 
63 61 17 58,2 
19 10 5 16,7 
Угол поворота коленвала при максимальном давлении в КС, 
град ПКВ после ВМТ ы ь о 
15 2 5 12,7 
59 73 29 57,4 
Максимальное давление, атм 69 77 31 66 
68,7 78,4 32,7 66,1 
ь 16 9 -13 12,4 
Угол поворота коленвала при максимальной скорости нарастания 
-13 2 -11 -11,3 
давления в КС, град ПКВ после ВМТ 10 -9 12 6 
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Окончание табл.2 









































1 2 3 4 5 

6,95 9,44 0,67 6,57 

Максимальная скорость нарастания давления в КС, атм/град ПКВ 11,8 9,8 8,5 11,3 
5,75 8,09 8,98 8,71 

м 25 10 23 23,3 

Угол поворота коленвала при максимальной средней температуре 23 6 18 208 
заряда, град ПКВ после ВМТ 15 2 22 14 4 
2240 2193 925 2103 

Максимальная средняя температура рабочего тела, К 2306 2175 942 2157 
2265 2135 1009 2126 

290 299 9 263 

Максимальная скорость нарастания температуры заряда, К/град 373 276 13 3273 
и 204 225 14 187 
0,05 0,03 0,23 0,09 

Доля теплоты цикла, выделившаяся до ВМТ 0,09 0,22 0,36 0,13 
0,12 0,59 0,36 0,19 

Доля теплоты цикла, выделившаяся до момента те В не р 
когда давление в КС максимально _. __ к — 
0,99 0,99 0,80 0,97 

Доля теплоты цикла, выделившаяся до момента, 099 008 0 76 098 
когда температура заряда максимальна 099 098 0.85 098 








Примечание: тактность двигателя - 4; число цилиндров - 4; рабочий объем камеры сгорания - 400 смз; частота вращения 
— 1000 / 2000 / 4000 мин”; степень сжатия - 10; отношение радиуса кривошипа к длине шатуна - 0,25; давление 
в начале такта сжатия - 1,0 атм; температура в начале такта сжатия - 300 К; начало впрыска топлива - 20 / 30 / 35 град 
до ВМТ; длительность впрыска - 35 / 30 / 40 град; момент поджига заряда - 7 / 13 / 9 град до ВМТ; средний диаметр ка- 
пель диз. топлива - 0,02 см; коэффициент диффузии паров диз. топлива - 0,1 см?/с; диаметр цилиндра - 8 см; ход порш- 
ня — 8 см; температура стенок камеры сгорания -— 600 / 650 / 700 К; средний коэффициент теплопередачи - 0,01 Вт/см?К; 
обратное время диффузии кислорода в факел - 1000 / 1700 / 2500 1/с. 


Приведенные в табл.2 данные согласуются с известными теоретическими положениями и экс- 
периментально установленными фактами. Развивае- 
мая мощность монотонно возрастает с частотой вра- 






















щения коленвала в интервале 1000- 4000 Р, ати Осредненное 
мин !, в то время как крутящий момент (он пропор- — 0,74 
ционален среднеиндикаторному давлению) максима- 75 — 10, 
лен в интервале 2000—4000 мин", Увеличение темпа ны 1,30 





вращения коленвала сопровождается возрастанием 
доли тепла, выделившегося до ВМТ с 9 до 19 %, 
что также отвечает современному пониманию ди- 
намики внутрицилиндровых процессов. 

Вклад трех макроскопических однородных 
зон заряда в силовую брутто-характеристику ДВС 
проиллюстрирован на индикаторной диаграмме 
(рис.2) для условий в средних строчках табл.? (но- 
минальный режим двигателя). 








Усовершенствованная модель удовлетвори- 
тельно воспроизводит кинетику СО-токсичности 
ДВС при изменении нагрузки, объема и закона по- 
дачи топлива. Кинетическая кривая [СО](Ф) для 
одноцилиндровой версии двигателя 10Д100 при 
увеличении мощности от соответствующей холо- 
стому ходу до номинальной показаны на рис.3. Раз- 
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Рис.2. Парциальные индикаторные диаграммы 
для различных зон ТВС двигателя ВАЗ-2108 
при частоте вращения коленвала 2000 мин" 
(при искусственно отключенном межзонном взаимо- 


действии) 
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гон ДВС в вычислительном эксперименте осуществлялся увеличением цикловой подачи топлива в 
4 раза; зависимость интенсивности впрыска топлива от угла поворота коленвала (ПКВ) не меня- 
лась. 


[СО], % 
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0 [2] Г) Г Т Г аа 
0 2000 4000 6000 8000 Ф, град. ПКВ 


Рис.3. Переходный процесс в одноцилиндровом двухтактном дизеле: 
кинетическая кривая [СО](Ф) 


Отображенный на рис.3 переходный процесс сопровождается увеличением удельного вы- 
броса СО примерно в 3 раза по сравнению с исходным и конечным установившимися режимами. 
Его продолжительность зависит от инерционности транспортного средства и динамических 
свойств трансмиссии. Разгон двигателя осуществляется за 20-30 оборотов коленвала, что соот- 
ветствует примерно двум секундам ускоренного движения городского автобуса. Важно отметить, 
что в городском ездовом цикле переходные режимы ДВС часто доминируют [1, 2]. 

Кроме того, реалистичные модели должны учитывать число цилиндров и тактность двига- 
теля, весь спектр выбрасываемых токсичных веществ (ТВ) и различия их индивидуальной ядови- 
тости. В результате усложняются кинетические кривые как отдельных ТВ, так и суммарной ток- 
сичности заряда [1, 2]. Отметим, что динамика интегральной токсичности ТВС и ОГ немонотонна. 

Рассчитанная согласно описанной модели, внутрицилиндровая кинетика углеводородов, 
альдегидов и МО, в т.ч. их содержание в ОГ тепловозного дизеля 10Д100 показаны на рис.4. 


-‹ [ВСНО], отн. ед. 
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Рис.4. Кинетика в среднеоборотных дизелях на примере тепловозного дизеля 10Д100: 
а- углеводородов; 6 - альдегидов; в — оксида азота 
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Условия расчетов соответствуют номинальной мощности для углеводородов и альдегидов, 
а для МО — номинальной (экстремальная кривая) и половинной (монотонная кривая) мощности. 
Содержание токсичных компонентов в ТВС дается в мольных долях. 

По завершении фазы активного горения основным углеводородом в ОГ является ацетилен 
(рис.4,а). Это объясняется наибольшей прочностью С=С связи ацетилена среди углерод- 
углеродных связей в моторных топливах. На этапе активного горения, однако, основным углево- 
дородом в заряде является этилен: его доминирование объясняется тем, что С!-углеводороды 
окисляются существенно быстрее, чем С›-улеводороды, а Сз-углеводороды, наоборот, заметно 
медленней. В результате С!-углеводороды сгорают практически мгновенно, а Сз-углеводороды 
в течение характерного времени окисления С›-углеводородов успевают разложиться термически 
на С!- и С›-фрагменты. 

Отметим, что рассмотренная ситуация отвечает условиям достаточно хорошего смесеоб- 
разования, при котором в выхлопе практически отсутствуют высшие углеводороды. Высокое ка- 
чество смесеобразования в рассмотренном примере выражается также в незначительном содер- 
жании альдегидов в выхлопе (рис.4,6)), что хорошо согласуется с современной теорией физико- 
химических процессов в ДВС [1, 2]. 

Экстремальный характер кинетической кривой № на рис.4,в, которая соответствует но- 
минальной мощности, означает, что при высокой теплонапряженности рабочего цикла ДВС кине- 
тика внутрицилиндрового оксида азота двухстадийна. На первой или горячей стадии концентра- 
ция МО успевает отслеживать среднюю температуру заряда. Вторая стадия реализуется, когда 
средняя температура становится слишком низкой для протекания реакций гибели МО, и его со- 
держание в ТВС «замораживается». 

При частичных нагрузках средняя температура рабочего тела в течение всего цикла не 
достаточна для «запуска» реакций гибели М№О, поэтому средняя концентрация окислов азота 
в рабочем теле монотонно увеличивается. В обоих случаях содержание МО в выхлопе близко к 
фактическому [1, 2]. Кроме того, рассчитанное соотношение МО-токсичности выхлопа при номи- 
нальной и половинной мощности получается вполне реалистичным. 

Еще один практически важный результат компьютерных экспериментов с моделями состо- 
ит в выявлении сравнительно компактного и универсального химического механизма горения 
нефтяного моторного топлива в современных поршневых ДВС, что исключительно важно для тех- 
нических приложений. 

Авторы, воспользовавшись разработанными ранее технологиями и средствами [10], вы- 
явили элементарные реакции, определяющие химизм рабочего процесса поршневых ДВС. Полу- 
ченная в результате оптимизированная схема превращений включает 70 реакций между 35 ком- 
понентами и пригодна как химический модуль перспективных моделей внутрицилиндровой дина- 
мики ДВС. 

Относительное отличие индикаторных (силовых, энергетических) характеристик, рассчи- 

тываемых по оптимизированной и полной схеме, не превышает 5 %. Отметим, что оптимизиро- 
ванная схема была получена путем объединения двух частично оптимизированных схем - для 
режима холостого хода и номинальной нагрузки. Существенно, что выделенные авторами ключе- 
вые реакции в большинстве своем признаны также определяющими горение нефтяного топлива в 
условиях жидкостных ракетных двигателей [11]. 
Выводы. Разработана динамическая модель внутрицилиндровых процессов поршневого ДВС, 
позволяющая рассчитывать как интегральные и мгновенные энергетические показатели, так и 
токсичность отработавших газов. Модель позволяет сконструировать минимальную кинетическую 
схему, адекватно отражающую динамические свойства современных поршневых двигателей. 
Указаны перспективные направления совершенствования разработки и возможности ее примене- 
ния для техники и технологий, сопряженных с использованием химических энергоносителей на 
транспорте. 
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